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ABSTRACT 

We describe how to probe new physics through examination of the form 
factors describing the Z—t—i couplings via the scattering process e~e + — > 
ti. We focus on experimental methods on how the top quark momentum can 
be determined and show how this can be applied to select polarized samples 
of ti pairs through the angular correlations in the final state leptons. We 
also study the dependence on the energy and luminosity of an e~e + collider 
to probe a CP violating asymmetry at the 10 -2 level. 
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1. Introduction 

Large CP violating effects are required to have the cosmological baryon asymmetry pro- 
duced at the weak phase transition!!] If such effects exist, can they be probed in hadron 
and electron colliders?! Schmidt and Peskini studied a CP violation mechanism suggested 
by Weinberg 4 ] and concluded that it is possible to probe a CP violating asymmetry in 
high-energy hadron collider experiments via the production of top quark-antiquark pairs. 
A detailed study on determining what observables are optimum for measuring the form 
factors governing the electric and magnetic dipole moments has been performed by At- 
wood and Soni.B In this paper it is shown that to probe the same level of CP violating 
asymmetry at the Next Linear Collider (NLC) with y/s = 500 GeV requires 3 x lO 4 ^ 1 to 
produce about 2 x 10 7 tt pairs, which is about a hundred times larger than the proposed 
integrated luminosity. (We use NLC to represent a generic e~e + supercollider B) To 
slightly improve the measurement on the CP violating asymmetry in e~e + collisions, we 
propose a different experimental observable from the one suggested in [Ref. 3| by measuring 



the polarization of the t and t in the process e _ e + — > tt — > bbl~l + vv. We show that the 
efficiency of untangling the different polarizations for the t and i increases as the energy 
of the e _ e + collider increases at the cost of decreasing the total event rate. The net effect 
is that a 500 GeV machine might do a better CP violation measurement than a 1 TeV 
machine. 

In present experimental data, the Z — t— i coupling was measured via its contribution 
to the decay of Z — > bb and the mass shift of the Z boson at the loop level. It was 
argued that based on the electroweak chiral lagrangian analysis, the form factors describing 
the Z — t—i interaction can be different from the ones predicted by the Standard Model 
(SM).@ We propose measuring these form factors by considering the polarization of the t 
and i. We show how to identify the polarization of the t and i and how to reconstruct 
the momenta of the t and t via the decay mode ti — > bl + vibl~vi, where / can be either an 
electron or a muon. We conclude that the form factors can be tested up to a few percent for 
the NLC with a luminosity of 50 fb _1 . For simplicity, in this paper we follow the argument 



in [Ref. 7| and assume that the dominant non-standard form factors are only present in 
the Z — t—i coupling and treat the couplings of the photon and boson with fermions 
as described in the SM. 

This paper is organized as follows. In Section 2 we present the interaction lagrangian 
and its corresponding helicity amplitudes which describe the production of ti pairs via 
e~e + annihilation using general form factors. There, we also present the Born level cross 
sections in the Standard Model. In Section 3 we review the requirements for discussing 
CP violation and in Section 4 we examine how well the form factors that govern the vector 
and axial vector interactions can be mneasured. In Section 5 we discuss how to reconstruct 
the top quark momentum in the e~e + — > ti — > bbl~l + uu and e _ e + — > ti — > I + jets mode 
including the effects of initial state radiation. In Section 6 we discuss top quark polarization 
and how to select polarized event samples, then in Section 7 we give our conclusions. 
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2. The General Form Factors and Helicity Amplitudes 



2.1. Conventions and Amplitudes 

In a recent paper! we studied the most general form factors for the coupling of t and i 
with either of the vector bosons 7 or Z. Decomposing this interaction over the basis given 

R L 

by the Clifford algebra in terms of the form factors i*\ 2 3 yields the general lagrangian, 



C'int —9 



W (f? {l) P- + f? (i?) p+) t - x -d v z^ v (F? {L) P- + pf (i?) P+) t 

+d"Z tM t (f^ L) P_ + F 3 Z(i?) P+) t + (F^ L) P- + if t (2.1) 

~d v A„ia^ (p 2 7(L) P_ + F 2 7(E) P+) t + d^A^i (p 3 7(L) P_ + F^ {R) P+ 



where P± = |(1 ±75), ia^ = -I[ 7 M y], an d v = (^Gp)- 1 / 2 ~ 246 GeV. 

Applying the Gordon decomposition, the vertex factors from the Z — t — t and ■y—t—i 
interactions become 

I* =|fly - BZ7 5 ) + PW "7 (iT (C, - B Z 7 5 ) + P(f) "; P(fT - F z ^]t, 
T ? =|f (r (A, - Brf) + "W-PW (C, - D 77 *) + + (B 7 - F 77 =)] ( , 



2.2 



where 



(2.3) 



. _ F h,Z](L) , F h,Z](R) 2m t ( F h,Z](L) , p [7,^](«) 
-^1 + *1 --^^2 +^2 

B b,z] - b \ ~ b \ 
r _ 2 ( F h,z](L) , F h,z](R) 

D M =- (F%ML)_ F b,z m 
E M J-2(F^ + F^), 
F [lZ ] =2 (pJ 7 ' Z](L) - F^ zm ^j . 

Since we ignore the masses of the incoming electron and positron, the P3 form factors 
in Eq. (2.1) (i.e., the P[ 7i z] an d ^Vy,Z] in Eq. (2.3), whose vector coefficient is the four- 
momentum of either the 7 or Z vector boson in our problem) do not contribute, however, 
we note that current conservation for the photon demands that 

Bi = z M^L ^m. F ^y (2 . 4) 
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In the SM, at the Born level, the form factors are 



A- 



3 



sin 9 



0, 



A' 



B z 



1 



2 cos 6w 
1 



sin 2 9 



w 



(2.5) 



2 cos 6*v(/ 
C Z = D Z = 0, 



where 0^/ is the weak mixing angle. Beyond the tree level, all of them except D, which 
controls the CP violation, have contributions due to loop corrections in the SM provided 
we ignore the small CP violating effects which reside in the Yukawa couplings that govern 
the interactions between the Higgs boson and the quarks. As argued in [Ref. 9| , the CP 
violating amplitudes in the Cabbibo-Kobayashi-Maskawa model are typically suppressed 
by a factor of order 10 -12 . 

The Z — e~ —e + and 7 — e~ — e + vertices are given by 



L Zee 



4p- + 4p+ 



and 



"/ee 



(e 7 L P- + elP+) , (2.6) 



where the SM values for these couplings are 
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1 



cos 9w 



sin 2 9w, and e 1 L = e 



R 



sin9 w . (2.7) 



The helicity amplitudes for e e 



tt are represented by (h e - , h e +, ht, %), 



where 



h e - = — , + respectively indicates a left-handed and a right-handed electron. Apart from 



the common factor 



2g 2 E 



with s = AE 2 , the nonvanishing helicity amplitudes from the diagram mediated by the 
Z-boson in e~e + — > tt are, in terms of these form factors Az, E>z-, Cz and DzP 

~-e L smd t [m t A z - K 2 C Z + EKD Z ] , 
= -e£(l + cos0 t ) [EAz + KBz], 
--e L {l-cosd t ) [EA z -KBz], 
~-e L sin 9 t [-mAz + K 2 C Z + EKD Z ] , 
~-e R sin 9 t [mtA z - K 2 C Z + EKD Z ] , (2.8) 
=e fl (l-cos0 t ) [EA Z + KB Z }, 
--- e R (l + cos 9 t ) [EAz-KBz], 
--e R sm9 t [-mtAz + K 2 C Z + EKD Z ] . 
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In the above formulas E = \^s/2 is half the center of mass energy of the e e + annihilation 
and K = \J E 2 — m\. The nonvanishing Born helicity amplitudes from the photon propa- 
gator diagram, (h e -, h e +, ht, hi)<y, can be easily obtained from the results in Eq. (2.8) by 
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setting e_c = e# = — sin 9yy and replacing the A, B, C, D form factors for the Z boson with 
those for the photon. The common factor in this case is 

2g 2 E 



The helicity amplitudes for the process e e + — > it are therefore obtained by summing the 
contributions from these two diagrams: 



(h e -, h e +, h t , hi) = 2g 2 E 



r(h e -,h e+ ,ht,ht) z , {h e -,h e +,htM) 



s-Ml 



+ 



(2.9) 



When calculating the cross section at the tree level, a color factor of 3 should be included 
for the top quark pair production. The spin average factor || should also be included for 
an unpolarized e~e + . 



Table 2.1: Born level production rates for e e + — > it at = 500 GeV 
and 1 TeV using m t = 140 GeV listed according to t and i polarizations. 



CROSS SECTIONS FOR VARIOUS t and t POLARIZATIONS 


tt Polarizations 


a(v^= 500 GeV) in fb 


a{y/s = ITeV) in fb 


RR = LL 


37.3 


2.69 


RL 


225 


66.7 


LR 


342 


104 


TT(TT) in = TT(U) in 


93.6 


29.8 


TT(|T) in 


290 


67.2 


2T(Ti)i„ 


164 


49.6 


TT(T|)± = TT(|t)± 


131 


32.5 


2T(TT)± = TT(U) ± 


190 


55.7 


t is R, i unpolarized 


262 


69.4 


t is L, i unpolarized 


379 


107 


top only is T(t)m 


258 


79.4 


top only is T(J.) in 


383 


97.0 


Unpolarized ti 


641 


176 
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2.2. Cross Sections in the Standard Model 

The production rates of polarized tt pairs at 500 GeV and 1 TeV e~e + colliders are 
tabulated in Table 2.T for a 140 GeV top quark. LL (RR) denotes the production rate 



of tit l ( tRtR ) events, where L labels a left-handed helicity and R labels a right-handed 
helicity. Likewise, RL (LR) denotes the production rate of tRti ( t^tR ) events. Besides the 
helicity cross sections, we have also included the cross sections for transverse polarizations 
of the t and i (indicated by T). In the first column of [table IO , the subscript in or 



_L respectively indicate whether the cross sections refer to transverse polarizations of the 
quark pair in the e~e + — > tt scatter plane or perpendicular to it. If the arrows are 
( ant i) aligned, the t and t spins are ( ant i) aligned. An upward (downward) pointing arrow 
in the case of polarizations perpendicular to the scatter plane means the spin is directed in 
the +y (— y) direction where y = z x p(£). (Recall the incoming electron beam is moving 
in the +z direction.) An upward (downward) pointing arrow in the case of polarizations 
in the scatter plane means the spin is directed perpendicular to the top momentum and 
toward the — z (+z ) direction. 

Since LR and RL only depend on the form factors A and B, it is possible to untangle 
the contributions from A and B by studying the polarized states tjJ,R and £r£z,. One 
important issue is how well the NLC can bound these form factors for a given integrated 
luminosity. 

If one can separate the polarization states of LL and RR from the others, then one 
can measure the asymmetry 

LL — RR 

A ' = llTrr (2 ' 10) 

or 

« LL — RR „ _ . 

A 2 = • (2.11) 

LL + RR + LR + RL K ' 

These asymmetries are important for studying CP violation. 



3. CP Violation 



3.1. Form Factors and CP Violation 



If the form factor D in Eq. (2.2) is not zero, the theory is CP violating!] Specifically, 
from Eq. (2.8) it is apparent that the difference between the RR and LL cross sections is 
linearly dependent on D through 



(RR - LL) oc Re 



A; 



s 



(3.1) 



where the "*" indicates complex conjugation and in this equation the photon is assumed 
to preserve its SM behavior, namely, _D 7 = 0. From the helicity amplitudes listed in 



LL-RR 



to measure 



Eq. (2.8), one can construct a quantity sensitive to CP violation, Ai — zl+rr j 
the real part of D.i The imaginary part of D can be examined by studying the transverse 
polarization of the top quark pairs perpendicular to the scatter plane.! The advantage 
of being able to measure the polarization states of the top quark pair is that one can 
measure Ai instead of Ai = LL+mi+^R+RL wmcn mav l° se some of its effectiveness due 
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Figure 3.1: This plot displays the variation of the asymmetry A\ in 
Eq. (2.10) with the mass of the Higgs boson for Weinberg's model using 
v/i = 500, 1000 GeV and m t = 140 GeV. 



Figure 3.2: This plot displays the variation of the asymmetry A\ in 
Eq. (2.10) with y/s for Weinberg's model using mt = 140 GeV and ran — 
0.1, 0.4, 1.0 TeV. 

to the dilution acquired by having a larger denominator from the RL and LR events when 
studying the CP violation effects. 

In Weinberg's model there is a dependence on the mass of the Higgs boson that enters 
in the form factor D, 

Re {B[7 , z]} = 2 I m{ Z 2 MB T (^y [l - Jl n (l + g) ], (3.2) 

where (3 = \J\ — Amf/s and A^^j 11 refers to the SM values at the Born level for A 1 and 
Az given in Eq. (2.5). For models with two or more Higgs doublets Weinberg showed that 
2Im{Z2} < \/2 with a reasonable choice of Higgs vacuum expectation values. When we 
compute A\ for a 500 GeV machine with m t = 140 GeV (see [Fig. 3.1| ) , we find A\ = —1.57% 



for mjj = 100 GeV and A\ = —0.43% for m# = 1 TeV. The percentage decrease in the 
magnitude of the CP violating asymmetry is less for a heavier top quark, but the general 
trend is a decrease in as the Higgs boson mass gets larger. 

Using Weinberg's model, it is apparent that there is also a variation of the CP violating 
asymmetry with the beam energy of the collider. In [Fig. 3.2| we plot Ai versus the center 



of mass energy for Higgs boson masses of 100, 400 and 1000 GeV using m t = 140 GeV. 
Comparing a 1 TeV machine with a 500 GeV machine we see that the A\ asymmetry 
value is slightly larger in magnitude for the higher energy machine. It is important to 
note, however, that the RR and LL event rates, which are the ones used for studying CP 
violation, are smaller by an order of magnitude for the 1 TeV machine when compared with 
the 500 GeV machine (see Table 2.1| ). Because of this reduction in statistics, the 500 GeV 



machine is preferred for testing CP violation. 

The question is, "What luminosity is needed to measure the CP violating asymmetry 
Ai at the 10~ 2 level at the NLC?" 
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3.2. Requirements for studying CP Violation 

To test CP violation, we use the dilepton mode, 

e~e+ -> it -> blri + uu, (3.3) 

which has a branching ratio of 2^ = ^- for Z = e, jx. At the Born Level in the Standard 
Model the unpolarized cross section for e~e + — > tt with rri( = 140 GeV is (see [Table 2.1| ) 



/_+ ^ J 641 fb for Js = 500 GeV, ,„ 
^ C 6 ^ tt ) = \l76fb for^i = lTeV, (3 ' 4) 

while the cross section for the production of t and t with identical helicities is 

f - + , n ( - + 4- J \ / 37 fb for v /i = 500 GeV, ,„ _x 

^(ee - fcfc) = ff (e e - = ( 2 .7fb for ^ = X Te V. " (3 ' 5) 

To measure an asymmetry A\ ~ 10~ 2 requires that the number of events, V, be large 
enough so that the statistical error in the measurement, -t= is less than 10~ 2 . For the 

optimistic value of ^= ~ 10 -2 we can estimate the luminosity required to observe this CP 
violating asymmetry, 

~ 2e [a (e~e+ -+ + ^ (e~e+ ^ t L t L )] ' (3 ' 6) 

where e is the efficiency of finding LL or ItLR events. Assuming we are able to solve for the 
kinematics of the dilepton mode in 70% of the events (see Section 5) and the acceptance of 
the kinematic cuts in selecting LL or RR events is 4 (see Section 6), then e ~ 0.18. This 
implies the required luminosity to observe CP violation through A\ ~ 10 -2 in top quarks 
of m t = 140 GeV is £ w 3 x 10 4 fb _1 at ^ = 500 GeV (i.e., about 2 x 10 7 it pairs are 
required) . 



4. How Well Can Form Factors be Measured at the NLC? 

The bounds determined from the experimental data, which limit the size of the form factors 
governing the top quark interaction with the vector boson, da not eliminate the possibility 
of nonuniversal couplings, as discussed by Peccei and Zhang,0 making it important to test 
the couplings as completely as we can. When determining form factors for an s-channel 
process where the intermediate state is a vector particle, one benefits bv the information 
conveyed by the variation in the angular distribution of the final state.E^I From Eq. (2.8) 
it is clear the it helicity states are produced with a polar angle distribution of either sin 2 6 
or (1 ±cos#) 2 , so noting that the form factors carry no angular dependences, it is possible 
to use the methods described above to solve for the direction of motion for the top quark 
and then use the differential cross section as described through the parameters co,c+,C- 
by 

da \/ E 2 — m'i 



± 

dcos6 32ttsE 



Co sin # + c+(l + cos#) z + c_(l-cos#n . (4.1 
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Figure 4.1: In da /d cos 9 we see the dependence of the production rates 
for polarized tt states on the top quark polar angle. This angle 6 is measured 
with respect to the incoming e~ beam which is taken to move along the +z 
direction. 



Figure 4.2: These plots present the angular correlation between the 
charged leptons produce from the polarized t and t decays in e~e 



+ 



tt 



bbl~l + uu. Boosting to the individual center of mass frame for the decays 
of the t and i, we plot the cosine of the angle between the momentum of 
each charge lepton and the boost axis as defined by the top momentum. 
Note how each combination of t,t helicities (RR,LL,RL,LR) occupy separate 
quadrants. 

From Eq. (2.8) and Eq. (4.1) we obtain 



co =co+ + Co-, 
c ± =6 (g 2 E) 



z m t A z - K 2 C Z ± EKD z 7 m t A 1 - K 2 C 1 ± EKD 1 



s — rrir 



+ 



z m t A z - K 2 C Z =f EKD z 7 m t A 1 - K 2 C^ q= EKD 1 
e R r, h e 



s — rrir 



J R 



(4.2) 



c± =3(g 2 E) 2 



Z EA Z ±KB Z J EA 1 ±KB 1 



s — rrir 



+ 



z EA z tKB z 7 EA 1 t KB 1 



s — rrir 



Untangling the form factors from such distributions is not trivial for various reasons. 
For instance, it is possible for there to be cancellations among the form factors making 
deviations from the SM difficult to observe. One result of these complications is that given 
a certain polar angle distribution described by a finite number of data points carrying some 
doubt in precision, there are different sets of values for the form factors that can satisfy 
the constraints such a distribution gives. With this in mind, nevertheless, we can ask how 
well experiments at the NLC might be able to constrain the form factors. 

We now proceed to ask how well the fit to the angular distribution can do to get the 
form factors. We show in Fig. 4.1| the lowest order SM production rates for polarized 
tt pairs as a function of the scattered polar angle, 9t, measured in the center-of-mass 
frame of the tt pair at the NLC. In Table 4.T| and [Table 4?2| we present bounds that 
represent a 68% and a 90% confidence level on the range for determiningthe form factors 
for the Z— t — t interaction. These bounds were obtained usina; MINUITEl to fit the polar 
angle distribution of the top quark as generated for the NLC (a/s=500 GeV, with 50 fb -1 
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Figure 4.3: The distribution for the production of polarized tt pairs, 
da/d[(l — r)/(l + r)] vs. (1 — r)/(l + r) for r = E(l + )/E(l~) when angular 
cuts are applied to the final state Z + and l~ in the respective decay frames 
of the t and t. 



Table 4.1: The upper and lower bounds indicate the statistical range 
within which we can determine the form factors from the angular distribu- 
tion of the top quark in e _ e + — > tt — > bbl~l + uu given a 50fb _1 luminosity 
at the NLC within a 68% confidence level. 



FORM FACTOR BOUNDS FOR THE NLC 


Form Factor 


SM VALUE 


Upper 


Lower 




0.395 


0.419 


0.370 


p Z(R) 


-0.175 


-0.153 


-0.197 


b 2 


0.0 


0.013 


-0.009 


b 2 


0.0 


0.013 


-0.009 



Table 4.2: The upper and lower bounds indicate the statistical range 
within which we can determine the form factors from the angular distribu- 
tion of the top quark in e _ e + — > tt — > bbl~l + vv given a 50 fb -1 luminosity 
at the NLC within a 90% confidence level. 



FORM FACTOR BOUNDS FOR THE NLC 


Form Factor 


SM VALUE 


Upper 


Lower 


r l 


0.395 


0.473 


0.311 


F Z{R) 


-0.175 


-0.107 


-0.254 


r 2 


0.0 


0.084 


-0.024 


F Z(R) 
b 2 


0.0 


0.084 


-0.024 



integrated luminosity) with no constraints on the kinematics. We expect 30,000 tt events, 
but if we focus on the mode where the two W's decay leptonically to ^ or there is a 
branching ratio reduction to 1,500 events. These 1,500 events were collected into twenty 
bins for the fit. Allowing only one form factor to vary at a time, the results presented in 
Fable 4.T| ( [Fable 4T2] ) indicate that within the 68% (90%) confidence limit, it should be 
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Figure 4.4: Same as [b ig. 4.1| except for a purely left-handed e beam. 



possible to find F± to within about 6% (20%), while 

F R 

can only be known to within 

roughly 11% (40%). The F2 form factors are zero in the SM, and the fit indicates that 
their values can be known to within about 0.0 ±0.01 (0.0^002) a ^ a confidence level of 68% 
(90%). 

By taking our knowledge of the polarization behavior to untangle some of the form 
factors before comparing them with data, it may be possible to improve these bounds. 
In particular from Eq. (2.8) it is apparent that the ti states of opposite helicity only 
depend on A and B. By making cuts (to be discussed in Section 6), we can isolate RL 
and LR contributions and investigate A and B from this data subset with a minimum 
contamination from the LL and RR states. Afterwards, we can proceed to the RR and 
LL contributions to study C and D. Performing the same procedures as we did with 
the unpolarized cross section, we take 750 events from the unconstrained LR sample and 
use MINUIT to fit for F^(Z) and find 0.395 ± 0.014 at a confidence level of 68%, which 
is roughly two percent better than we did with the unpolarized distribution. Further 
improvement in the form factor determination can come from the increased statistics (by 
about a factor of 7) obtained by including the events from e~e + — > ti — > I + jets in the 
analysis. For the I + jets mode, where the branching ratios take 9, 000 unpolarized events 
from the original 30, 000, the results indicate that within the 68% (90%) confidence limit, 

it should be possible to find f[ L ^ to within about 3% (8%), while F± can be known to 
within roughly 5% (18%). In this case the F2 form factors can be known to within about 
0.0 ± 0.004 (0.0 ± 0.02) at a confidence level of 68% (90%). 

More information can be gained by focusing on polarization states for the top quarks 
besides pure helicity states. For example, the interference terms that result between the 
product of the imaginary parts of the amplitudes with the real parts do not exist in the 
total rate expressed by Eq. (4.1) and Eq. (4.2), but they do appear in the rate for top quark 
spins perpendicular to the scatter plane. To see this, note that the polarization asymmetry 
comparing top spins pointing in opposite directions perpendicular to the scatter plane is 
proportional to 



J2 Im 



(h e ^,h e+ ,h t = +, hi) (h e -,h e+ ,h t = ~,h^ • (4.3) 



Thus, angular correlations and decay plane correlations will further improve the determi- 
nation of the form factors F\ and F%. 

Even more useful, however, may be the purity of the ti polarization which is obtained 
by polarizing the electron beam.113 The majority of the cross section is comprised of 
left-handed top quarks from the start. By selecting a left-handed polarization for the 
electron beam, helicity conservation suppresses the production of right-handed top quarks 



12 



that move along the initial electron direction of motion. This can be seen in [Fig. 4.4 
where for each top quark helicity combination we plot da /d cos 6 vs. cos#. What this 
provides is a cleaner measure of the form factor dependence EA + KB, as apparent from 
the amplitudes shown in Eq. (2.8). Similarly, a right-handed polarization for the electron 
beam will provide a cleaner sample of right-handed top quarks where the form factor 
probed is EA - KB. 



5. How to reconstruct the Kinematics of the Top Quarks 

In reality, the initial state electron or positron at high energies will radiate photons along 
the beam axis either due to initial state radiation (ISR) during the e + e~ interaction or due 
to the interaction with the classical electromagnetic field formed by having a dense bunch of 
charged particles in the interaction region (beamstrahlung) . Such radiation tends to move 
along the beam direction such that neither the center-of-mass energy nor the boost of 
the ti pair is known, complicating the analysis discussed in the previous section. Despite 
this, the effects of initial state radiation are such that most of the time the center-of- 
mass energy of the ti pair is close to the beam energy. Furthermore, it is possible to 
design experiments such that the beamstrahlung is minimized, so we only consider the 
bremsstrahlung effect. Nevertheless, it is still possible to determine the momentum of the 
top quark. 

We direct our attention to solving for the top quark momentum for two cases: The 
first case will be the mode where the top quark decays leptonically while the top antiquark 
decays hadronically, e~e + — > ti — > bl + ubqq'; and the second case will be the mode where 
both top quarks decay leptonically, e~e + — > ti —>■ bl + vbl~v. For the initial state radiation 
we use the distribution for gettingan e~ (or e + ) with a fraction z of its original beam 
energy, given by Kuraev and Fading to order A, 

D (z, s) = ~(l-z) (l + |f) - \ X (1 + *) > 

where A = ^(ln^j — 1) for original center of mass energies of y/s and electron mass m e . 
In Section 6.2 we show that there was little difficulty in associating the proper bottom 
flavored quark with the proper top flavored quark, so we will drop the concern about this 
ambiguity for this part of the discussion. 

5.1. e~e+ -> ti -> bl+vbqq' 

When one top quark decays hadronically and the other decays leptonically, the only infor- 
mation we are missing to describe the momenta and energies of the event completely is the 
component of the neutrino momentum that moves along the beam axis. The transverse 
momentum of the neutrino is known by applying conservation of momentum. Next, we 
apply a familiar formula which provides a doubly degenerate solution for the longitudinal 
momentum of the neutrino, 

a ± ,/a 2 - (1 - b 2 ) (a 2 -b 2 p 2 (u)) 

P* iy) = — v - ^rp • 
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Figure 5.1: The center-of-mass energy of the e e + annihilation when 
ISR is taken into account. 

In the above formula, 

-(p r (6/+)-p r (,) + iK-Mg + )) E(U+) 

a ~ pAbi+) ' b -^wy (5 - 2) 

where p(z') = (p T (u),p z (u)) is the momentum of the neutrino separated into directions 
respectively transverse and along the Z-axis, while p(bl + ) = (E(bl + ), p(bl + )) describes the 
four-momentum of the two-particle system of mass p(bl + ) -p(bl + ) = M^ [+ composed of the 
bottom quark and the charged lepton. With the idea that polarization studies are a means 
of investigating top quark properties after its discovery, the top quark mass is assumed to be 
known and its decay width narrow. In our Monte Carlo study, the widths of the W bosons 
have been included for event simulation, so this computation was performed by scanning 
a range of values for m(W~) and m(W + ) about the peak value of m\y and accepting 
the first solution that was found. Attempting to account for some detector inefficiency, 
this procedure also included a gaussian smearing of the momentum for the b and b with 
AE(b)/E(b) = 0.b/y/E(b) and for the visible leptons with A J E(/ ± )/E(/ ± ) = 0.15/ ^WW). 
(More detailed discussion is given in Section 6.2.) Random chance yields a 50% probability 
for selecting the correct solution for p z {y\ but we need to do better than that. What needs 
to be resolved then, is how to remove this degeneracy in the ^-momentum of the neutrino. 

Table 5.1: The efficiency in solving for the top quark momentum us- 
ing initial state radiation, smearing, and a Breit-Wigner width on the 
VF-bosons is described here by three values: the percentage of solveable 
configurations; the cosine of the angle between the true top quark direction 
and the direction given by solving the kinematics (£ = cos(^ rue , t so i ve d))', 
the root mean squared of (. 



SOLUTION EFFICIENCY(unknown m w , smearing and ISR included) 


l + + jets 


%Solved 


<c> 


< c 2 > - < c> 2 


Best mw 


84 


0.978 


0.090 


Best E cm 


72 


0.981 


0.073 


l+l-bbvv 


%Solved 


<c> 


< c 2 > - < c> 2 


Best E cm 


68 


0.957 


0.142 



We have two straightforward means at our disposal for selecting the best neutrino 
solution. In one case we take advantage of the fact that the width of the TU-boson is 
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narrow with respect to its mass and select the p z {y) that gives a value for (p(u) + p(e + )) 2 
that is closest to rri^r. In the second case we take advantage of the feature that the 
initial state radiation is strongly peaked for soft emissions, as shown in |Fig. 5.1| , and select 
the Pz(y) that gives a value for the center-of-mass energy that is closest to the original 
beam energy. To measure how these two methods compare for determining the top quark 
momentum, we tabulate three values in |Table 5.1|. First, as a matter of statistics, we give 



the percentage of the time this method was able to find a solution. Selecting the best m\y 
provided the largest number of solutions where only 16% of the cross section was lost while 
choosing the smallest loss of energy through radiation failed to find a solution 28% of the 
time. What may be more important in polarization studies than losing an extra 12% of 
the events, however, is getting a more precise measure for the direction of the top quark. 
Averaging the cosine of the angular separation between the true top quark direction and 
the direction obtained by selecting the solution to the quadratic formula for p z {y)-, which 
we label C = cos(ttrue,t solved), it was found that the minimum radiation method had a 
value for ( slightly closer to unity. Not only is the average of this angle important, but 
also the spread of this distribution, < ( 2 > — < £ > 2 , where it was better to accept those 
solutions which provided the value of M(tt) closest to the original beam energy. 

With either of these two methods, there is little loss in statistics and very little change 
in the polar angle distribution of the top quarks, making it possible to study the form 
factors through the polar angle distribution of the top quark. 



5.2. e-e+ -»• tt -> bl+vbl 



v 



When the top quark and antiquark decay leptonically and both the boost and the mass of 
the tt pair are unknown, there is a great deal of information that must be extracted out of 
the measurement of the four visible particles (6, b, Z + ,Z~ for I = e,/i). Detailed discussion 
on such kinematics can be found in Kef. 16| . The four-momenta of both the v and v must 



be determined using the center of mass energy as an unknown. This means we have eight 
unknowns. Generalizing the procedures that have been discussed here ealier, however, we 
find that we can solve this system exactly, albeit within an eight-fold degeneracy at most. 

Determining the neutrinos' momenta is accomplished by utilizing the fact that this 
process contains the decay of four particles whose masses either are or will be known and 
whose decay widths are narrow. From the start, we have the conservation of transverse 
momentum, 
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v T + pt iy) + PT (v) = 0, (5.3) 



where we define the four-vector for the sum of the energies and momenta of the visible 
particles to be 

V = (V (0) , V) =p(b)+p(b)+p (e+) + p (e~) . (5.4) 

The conservation of transverse momentum provides two equations in pursuit of solving for 
the p(v) and pip)- Next, from the the decays of the W's and t's and the masses of the v 
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Figure 5.2: The boost, da /d(3 vs. (3 = ^fhjnjgm , as generated by ISR. 



(5.5) 



Figure 5.3: The angular distribution da /d cos 9 vs. cosO as dictated by 
the top momentum provided by the Monte Carlo compared against the top 
momentum found when solving the kinemtics incuding the effects of ISR. 

and v we have six more relations, two of which are nonlinear: 

2 2 

Top quark mass: m (t) = (p (b) + p(v) + p , 

2 — 2 

Top antiquark mass: m (t) = (p (b) +p(v) +p (e - )) , 

2 2 

W + mass: m (W + ) = (p (e + ) + p (z/)) , 

2 2 

W~ mass: m (W~) = (j> (e _ ) + p (u)) , 

2 — * — * 

neutrino mass: = E (u) — p (z/) ■ p [v] , 

2 -* — » 

antineutrino mass: = E [p) — p (u) ■ p (p) . 

In this study we have used the narrow width approximation for the top quark and top 
antiquark. 

These relations provide two coupled quadratic realtions in E{v) and E(u). Combining 
these two quadratic equations provides a quartic equation in E(u) which can be solved 
exactly, yielding a four-fold degeneracy at most for the neutrino energy. Substituting each 
solution for the neutrino energy, one at a time, into either of the two quadratic relations in 
E(u) and E{y) provides two solutions for E{v) for each E(u). The result is at most eight 
pairs of solutions for (E(u),E(u)). Since the remaining equations are linear, the solutions 
for the momentum components are unique for each given (E(u),E(u)) pair. What remains 
is to decide which solution is the correct one. 

Since the initial state radiation is peaked for soft emissions, as shown in [b'ig. 5.1| , the 



method we use for choosing among the degenerate solutions in this top decay mode is to 
select the set of that provides a value of M(tt) closest to the original beam 

energy. In [1'able 5.T] , we see that we were able to obtain solutions for 68% of the cross 



section with < C >= 0.957 and < C 2 > - < C > 2 = 0.142. 
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Figure 6.2: E(l + ) and E(l ) distributions for e e + — > tt when both final 
state top quarks have identical helicities. 



Figure 6.3: The distribution for the production of polarized tt pairs, 
da/d[(l - r)/(l + r)} vs. (1 - r)/(l + r) for r = E{l + )/E{r). 

5.3. Initial State Radiation 

The effects of ISR on the angular distribution of the top quarks is minimal. To demonstrate 



this, we show in |Fig. 5.1| and [Fig. 5.2| that the initial state radiation is soft, keeping the 
energy of the process quite close to the beam energy for a majority of the events. The 
tt system is hardly boosted as illustrated in Fig. 5.2 , where (3 = Wrr^rx ■ When we 



E(t)+E(t)- 

take events generated by our Monte Carlo and compare the angular distributions between 
the kinematics obtained by solving the eight simultaneous equations and the actual value 
generated by the program, we find only a slight difference, as shown in [Fig. 5.3| . 

With the t and i momenta determined, it becomes feasible to select a sample of tt pairs 
that is strongly biased to a particular polarization state. 



6. Determining the Polarizations of the t and the t 



6.1. Polarized tt Production 

Here we consider the dilepton mode. We showed in [Ref. 8| that the polarization of the t 



(and the t) can be self-analyzed from the decay t — > bW + — > bl + ui (i — > bW~ — > bl~i>{). 
For a heavy top quark, the preferred moving direction of the l + and the l~ in the center-of- 
mass frame of the tt pair are shown in Fig! 6.1| . Because of the correlations, it is possible to 



distinguish different polarimtion states of the tt pairs by the energy distribution, da /dE{l), 
of the leptons l + and Z~.Q'EJ'EJ Because the asymmetry in E(l + ) and E(l~) is magnified 
when the t and the i are boosted, we anticipate that it becomes easier to distinguish 
different polarization states of the tt pair when the energy of the e~e + collider increases, 
but it must be kept in mind that the production rate of tt decreases with higher beam 
energies. At the NLC, the energy distributions of the l + and l~ are shown in |Fig. 6.2| . 



As discussed in the previous section, Ai is cleaner than A2 because it does not contain 
the relatively large contributions from LR and RL in its denominators. Obviously, this 
can be useful only if one can separate the various polarization states for the tt pairs. One 
method to do this which applies the physics contained in the energy distributions of [Fig. 6^2 



is to make a selection of events using the energy ratio r = E(l + ) / E{1 ). The events with 
r>l are mainly RR events, r<l are mainly LL, and r ~ 1 characterizes both RL and 
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LR events. In order not to induce the CP asymmetry while applying the kinematic cut on 
this variable r, it is necessary to make a symmetric cut on r. For instance, one can select 
the RR events by requiring r > 3 and the LL events by requiring r < 1/3. In [Fig. 6.3] , we 
show the distribution of (1 — r)/(l + r) for various polarization states of the ti pair. 

Another method we propose is based on the fact that the moving direction of the l + 
and l~ are strongly correlated in the center-of-mass frame of the ti pair. To illustrate this 
point, we show in |Fig. 472 the strong correlation between cos#^+ and cos 0^-, where t [± is 
the angle between boost axis of the t and and the momentum of the l + (Z _ ) in the rest frame 
of the t (i). The four populated corners correspond to the four different helicity states of 
the ti pair. Therefore, one can select individual polarization states of ti by making cuts 
on this plot. However, it is important to make sure that the kinematic cuts won't induce 
an artificial CP asymmetry. For example, one can select the RR and LL data sample 
by making a square cut of equal area on the cos6>^+ > 0, cos0 ti - > and cos6>^+ < 0, 
cos0 t i~ < corners without inducing an artificial CP asymmetry when studying A\. The 
cos t i+ > 0, cos t i~ < and cos ti + < 0, cos ti - > corners are suitable to measure 
the form factors A and B from the RL and LR data sample. Similarly, information on 
the form factors C and D are contained in the LL and RR samples. The key element of 
this approach is to be able to determine the moving direction of the t via the decay mode 
of ti — > bl + V[bl~vi. This method is therefore useful even when some kinematic cuts are 
being imposed on the final observable partons as usually done in reality for detection and 
triggering. In [Fig. 4.3| we show the effect of selecting various corners in the configuration 
space of Fig. 472] on the distribution in (1 — r)/(l+r). The acceptance of the kinematic cuts 
in selecting LL or RR events is about \. Different polarization states of the ti dominate in 
different regions of (1 — r)/(l + r). Therefore, it is possible to enhance the probability of a 
particular polarization of the ti pair by making cuts on (1 — r)/(l+r) as to be performed 
in the following analysis. We also note the relative enhancement of the RR and LL rates 
displayed in Fig. 4.3| result of these selections. 



6.2. The Bottom Quark Ambiguity 

If the b is not distinguished from the b, there are two solutions provided by this procedure 
because it becomes possible to assign mistakenly the b with the t. Nevertheless, this causes 
only minor difficulties. To break the degeneracy, one first solves the kinematics given one 
combination of jet momenta assignments and then repeats the procedure for the other 
combination, then one compares the two results. 

We performed a simplified computation without ISR where it was not necessary to 
apply both of the constraints provided by m(W + ) and m(W~), since then we know the 
boost of the system. Furthermore, with the center-of-mass energy known, we did not need 
to invoke the nonlinear relations restricting the neutrinos to be massless. This means that 
given one W boson mass, we were able to solve for the other. The tendency is that the 
combination where the b quark was misidentified as the b quark generates unreasonable 
solutions for the W boson mass that was left as a free parameter, allowing a proper 
identification of the b quark jets. In particular, it was chosen to drop any solutions that 
gave a negative W boson mass (this condition may be tightened) . If both solutions yielded 
positive masses for the gauge bosons, they were subjected to a further selection criteria 
regarding the neutrino masses. 
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One effect of using a value for the mass of the W~ different than ■sj (p(e~) + p(v)) 2 in 
the above procedure is that the solution produces magnitudes for the neutrino momenta 
which are no longer equivalent to their own energies. To remove some of the inefficiency 
caused by the finite width, T\y, two solutions were determined. The first solution was 
obtained by fixing the mass of the W~ boson to m\y and solving the kinematics. For the 
second solution we left the mass of the W~ boson free and fixed the mass of the W + boson 
to rri]y. An artifiact of this procedure is that the worst solution of the two tends to provide 
neutrinos with larger values of p(^) 2 and p(v) 2 - Since the neutrino is massless, selecting the 
solution that gives the neutrino masses closest to zero ameliorates the problems introduced 
by r^/ and the remaining problem introduced by the inability to distinguish the b quark 
from the b antiquark. 

We have used a Monte Carlo program to test the success rate of this algorithm and 
found a 99% efficiency in reconstructing the moving direction of the t by assuming a hun- 
dred percent detection efficiency of the energy of the b and b. About 0.5% of the solutions 
misidentified b quarks as b antiquarks and the other 0.5% consisted of kinematic configu- 
rations which could not provide a solution given our criteria. To determine whether one 
can measure CP violation effects better than a percent from this method requires further 
study to ensure that no artificial CP violation is induced by removing from consideration 
those events for which no solution was found. Attempting to account for some detector 
inefficiency, this procedure was also studied with the inclusion of a gaussian smearing of 
the momentum for the b and b with AE(b)/E(b) = 0.5/ a/ E(b) and for the visible lep- 
tons with AE(l±)/E(l±) = 0.15/ y/E(l±). The results in this case gave a 96% efficiency 
in assigning the bottom momenta to the proper top quarks where 3.3% of the time this 
assignment was incorrect and 0.7% of the events provided no solution given our criteria. 
It is expected that the efficiency in removing the bottom quark ambituity should improve 
slightly by selecting the best values of piy) 2 and p(v) 2 among a range of choices where we 
do not restrict ourselves to fixing the W boson masses to m^y, but rather allow ourselves 
to check a variety of solutions for W boson masses within, e.g., mw ± ^w- Note that 
in distinguishing the b from the b through their kinematics as we have described, we have 
not taken advantage of other information whichcan be made available such as the charge 
of the leading particle in the bottom quark jet.13-3 

Summarizing, we described how the top quark analyzes its own polarization through 
its decay t — > bW + — > b l + u such that a moving top (anti) quark with right-handed (left- 
handed) helicity tends to produce the Z + (Z - ) with its momentum aligned with the top 
(anti) quark direction of motion. We have demonstrated that the top quark momentum 
can be determined well including smearing and width effects with very little ambiguity in 
deciding which bottom flavored quark goes with which top flavored quark. If the CP viola- 
tion effects are very small, however, the error introduced by the ambiguity in distinguishing 
the b from the b may interfere with the measurement. 

A 1 TeV machine can do better than a 500 GeV machine regarding the determination 
of the top quark momentum, polarization, and which bottom quark should be associated 
with the t, all due to the larger boost the top quark receives. 



19 



7. Conclusion 

We have shown how to determine the top quark momentum in e~e + — > tt — > bbl~l + vv and 
e~e + — > tt — > / + jets by solving for the momenta of the neutrinos, including the effects 
from initial state radiation. Assuming the mass of the top quark is known, an accurate 
determination of the top quark moving direction was obtained, even when including de- 
tector effects by smearing the particles' momenta and taking into account the finite width 
effects from the W boson. The efficiency of this algorithm was affected very little by the 
degeneracy created by treating the b and b as indistinguishable, and it is possible to op- 
timize the efficiency to account for width effects in the mass of the W boson. With the 
top quark direction known, the door is open to performing studies using the polarization 
of the top quark. 

Through angular correlations between the lepton and top quark directions, we have 
shown how to separate the different polarization states for the tt pair. Applying these 
tools, we have examined how CP violation effects and form factor magnitudes may be 
studied in e~e + — > tt — > bbl~l + vv and concluded that large luminosities (3 x lC^fb" 1 ) are 
required to measure a CP violation asymmetry Ai at the order of 10~ 2 at the NLC. On 
the other hand, the form factors can be measured to some accuracy. For the + jets mode 
at a 500 GeV machine , the results indicate that within the 68% (90%) confidence limit, 

it should be possible to find f[ L ^ to within about 3% (8%), while f[ R ^ can be known to 
within roughly 5% (18%). In this case the F2 form factors can be known to within about 
0.0 ± 0.004 (0.0 ± 0.02) at a confidence level of 68% (90%). Isolating the RL and LR 
contributions can improve the precision of the measurements of these form factors. For 

instance, we saw that in the dilepton mode that can be measured to about 6% at the 
68% confidence level using unpolarized beams, yet by selecting the RL contribution using 
the methods discussed in this paper, it appears possible to improve this bound to 4% in 
this mode. 

A 1 TeV machine can do better than a 500 GeV machine in determining F 1 ' because 
the relative sizes of the RR and LL production rates become small and furthermore, the 
top quark is boosted more in a 1 TeV machine thereby allowing a better determination 
of its direction. A 1 TeV machine makes the CP asymmetry measurement more difficult, 
however, because of the smaller rates in RR and LL. 
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Table Captions 

Table 2.1 Born level production rates for e~e + — > tt at y/s = 500 GeV and 1 TeV using mt = 
140 GeV listed according to t and i polarizations. 

Table 4.1 The upper and lower bounds indicate the statistical range within which we can deter- 
mine the form factors from the angular distribution of the top quark in e~e + — > tt — > 
bbl~l + vv given a 50 fb -1 luminosity at the NLC within a 68% confidence level. 

Table 4.2 The upper and lower bounds indicate the statistical range within which we can deter- 
mine the form factors from the angular distribution of the top quark in e~e + — > tt — > 
bbl~l + vv given a 50 fb^ 1 luminosity at the NLC within a 90% confidence level. 

Table 5.1 The efficiency in solving for the top quark momentum using initial state radiation, 
smearing, and a Breit-Wigner width on the VF-bosons is described here by three 
values: the percentage of solveable configurations; the cosine of the angle between 
the true top quark direction and the direction given by solving the kinematics = 
cos(ttrue,t solved))', the root mean squared of (. 
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Figure Captions 

Figure 3.1 This plot displays the variation of the asymmetry A\ in Eq. (2.10) with the mass of 
the Higgs boson for Weinberg's model using y/s = 500, 1000 GeV and mt = 140 GeV. 

Figure 3.2 This plot displays the variation of the asymmetry A\ in Eq. (2.10) with v / s for Wein- 
berg's model using mt = 140 GeV and m# = 0.1, 0.4, 1.0 TeV. 

Figure 4.1 In do jd cos 9 we see the dependence of the production rates for polarized tt states on 
the top quark polar angle. This angle 9 is measured with respect to the incoming e~ 
beam which is taken to move along the +z direction. 



Figure 4.2 Same as |Fig. 4.1| except for a purely left-handed e beam. 

Figure 5.1 The center-of-mass energy of the e~e + annihilation when ISR is taken into account. 
Figure 5.2 The boost, da/d/3 vs. (3 = ^m^^m > as generated by ISR. 

Figure 5.3 The angular distribution da jd cos 9 vs. cos9 as dictated by the top momentum provided 
by the Monte Carlo compared against the top momentum found when solving the 
kinemtics incuding the effects of ISR. 

Figure 6.1 Helicity arguments control the preferred moving direction for the e + ,e~ as they are 
produced in the respective rest frames of the t, i decays (t — > be + v, t — > be~&). The 
subscript R (L) indicates right-handed (left-handed) quark helicity. 

Figure 6.2 E(l + ) and E(l~) distributions for e~e + — > tt when both final state top quarks have 
identical helicities. 

Figure 6.3 The distribution for the production of polarized tt pairs, da/d[(l — r)/(l + r)] vs. 
(1 -r)/(l + r) for r = E{l + )/E{l~). 

Figure 6.4 These plots present the angular correlation between the charged leptons produce from 
the polarized t and t decays in e~e + — > tt — ► bbl~l + vv. Boosting to the individual 
center of mass frame for the decays of the t and i, we plot the cosine of the angle 
between the momentum of each charge lepton and the boost axis as defined by the 
top momentum. Note how each combination of t,i helicities (RR,LL,RL,LR) occupy 
separate quadrants. 

Figure 6.5 The distribution for the production of polarized tt pairs, dcr/d[(l — r)/(l + r)] vs. 

(1 — r)/(l + r) for r = E(l + )/E(l~) when angular cuts are applied to the final state 
l + and Z - in the respective decay frames of the t and t. 



